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FOREWORD 


The subject of NASA Research Grant NGL 05 - 010-176 is the theoretical 
and experimental study of space-related plasma wave propagation and 
resonance phenomena. Research supported by the grant has been proceeding 
under the direction of Prof. F. W. Crawford since its Inception on 
1 December 1966 . This is the sixteenth semiannual report, and covers 
the period from 1 July to 31 December 197 1 *. 
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I NTROLlt’CT I ON 



The research to be discussed In this semiannual report Is concerned 
primarily with plusmu wave nnd resonance phenomena related to the Iono- 
sphere. Much of our previous work under the grant has Involved direct 
laboratory simulation of puzzling effects occurring In space plasmas. 

Such projects led to better understanding of "Alouette" resonances, the 
resonance probe, some aspects of whistler propagation and triggered VLF 
emissions, and contributed to the development of a variety of new plasma 
diagnostic techniques. Particularly during the last decade or so, know- 
ledge of plasma wave and resonance phenomena, In both the linear and 
nonlinear regimes, has advanced rapidly as a result of considerable 
study worldwide of space and laboratory plasmas, and many of the puzzling 
phenomena have been explained. It is consequent ly possible to plan 
experiments with greatly increased confidence In the theoretical pre- 
dictions. The motivation ard aims of our research program have progressed 
with these advances, and bring us naturally to the projects to be 
discussed in Sections II and III, particularly in relation to the NASA 
Space Shuttle Program. 

The Space Shuttle Program offers exciting new possibilities for 
plasma physics experimentation. It is expected that, beginning in about 
I98O, the Shuttle Orbiter will carry as one of its regular payloads a 
"Spacelab", to be constructed by the European Space Research Organization. 
The Principal Investigator for this grant is a member of the Atmospheric, 
Magnetospherlc, and Plasmas-in-Space (AMPS), Science Working Group, set 
’•o by NASA to study the appropriate instrumentation for Spacelab, and 
'.cts as chairman oi the Plasma Wave Phenomena Section. 

If the AMPS enterprise is to be successful, new space plasma experi- 
ments should be submitted to critical scrutiny for scientific merit and 
experimental feasibility. To be attractive, they should profit from the 
unique plasma environment and parameter ranges accessible to the Space 
Shuttle to perform basic plasma experiments not feasible in the labora- 
tory, or to greatly improve others by virtue of the large volumes of 
essentially fully-ionized, c >llisionless plasma available for study; 
they should provide refined diagnostics for plasma conditions local to 
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the vehicle, und they should develop new techniques for Improving under- 
standing of the ioi '^sphere itself by strong perturbation and by nonlocal 
probing, e.g. with electron beams or waves. Our alms under this grant 
are to test and refine various ideas for space plasma experimentation 
using Spacelab. Over the last twelve months, three Ph.D theses have 
been completed with support from the grant. This has given us the 
opportunity to phase out a number of projects, and to initiate others of 
more direct relevance to the Shuttle. The current program is described 
in Section II. Brief comments on likely directions for the second half 
of the grant year are made in Section III. A bibliography of Reports, 
Conference Papers, and Publications resulting from the grunt during the 
reporting period, 1 July - 31 December 197*S * 8 Kiven in Section IV. 
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II. CUKKKNT RESEARCH PROGRAM 

A. Whistlers and Alfven Woven 

We have been studying various aspects of whistler propagation over a 
period of several years, beginning with laboratory measurements of 
whistler dispersion characteristics in a cold colllslonul plasma, con- 
tinuing with numerical studies of small-slgnnl whistler Instabilities, 
for propagatlo. both parallel' ' and oblique to the earth's 

magnetic field In nonthermul plasmas, and finally extending this work 

21 22 

to the nonlinear phenomena associated with triggered VLF emissions, * 

2 3 

modulatlonal instabilities, and Alfv 6 n wave excitation. We have also 

10 

studied excitation of VLF waves by charged particle beams. * and ULF 

24 

waves by metallic electric and magnetic dipole anten* s, as >art of 
our most recent work. During the reporting period, except for some minor 
revisions to a previously published paper, and to a manuscript submitted 

for publication 1r> the J. Geophys . Res . during the last reporting 

2 1 * 

period, the emphasis has been on whistler excitation by electron and 

. , ; » 10 

proton beams. 

Beam Excitation of Whistlers : A new method of generating VLF signals 

has been proposed recently In which helical electron beams, fired from 

OK 

satellites, are to be used to radiate whistlers . 1 The satellites are 
assumed to bo instantaneously on the L - 4 shell at an altitude of 
500 km. Electron beams of about 1 A current and 10 keV energy are to be 
Injected from 20 satellite-mounted guns forming a circular array of the 
order of 100 m diameter. The frequencies and amplitudes of the whistlers 
radiated will depend on the pulse frequency, and angle at which the 
beams are injected Into the earth's magnetic field. It has been claimed 
that the total power radiated from such a set-up will be at least of the 

-4 2 

order of 10 W, and that the power flux will be of the order of 10 W/m‘ . 
The former is comparable to that achievable by large ground-based VLF 
transmitters. 


* References I-I 5 arc to be found in Section IV. The remainder are 
given in Section V. 
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W r have been considering the proposul carefully, and trying; to 
determine the feasibility of currying it out with Spacelab. We noted 

14 

in SAR r) that the proposer of the idea had neglected to consider the 
stability of the beam to longitudinal wave growth. According to rough 
calculations of the growth rate of longitudinal instabilities, due to 
beam-plasma interaction, we believe that the beam cannot travel more 
than nbout 10 m before its helical s ructure and phase coherence are 
broken up by such instabilities. Thi \ would have the effect of reducing 
the powc r radiated to about four crders of magnitude below the 10 W value 

mentioned above. The power calculated by us was that radiated within a 

o 

propagation-vector cone-angle of JO , at fixed frequencies of 1C and 

4^ kHz. 

During the reporting period, we have refined our work by calculating 
the power radiated for a wider range of f requencies, and from within 

propagation-vector cone-angles corresponding to a range of duct enhance- 

27 

ments (due to plasma density inhomogeneity' ) between 1.02 and 1.2‘ . 

These more refined results verify that the radiated power is expected 
to be of the order of several mW. In addition to beam-plasma interaction, 
we have considered the effects of diffusion, space-charge spreading, and 
collisions, on the power radiated. We have also considered the merits 
of using an ion beam in place of the electron beam, and found the former 
to be superior to the latter with regard to radiation efficiency, 
although at the cost of a somewhat increased accelerating potential. 

We have identified in our analysis the term which gives ri *e to wave- 
particle interaction. This has allowed us to show that the waves which 
undergo wave-particle interaction are well outside the range which can 
be ducted, and therefore to explain why the power calculated by us is 
so much smaller than that calculated by Dowden. The results obtained 

have been described in detail in a recent report which has been 

10 

submitted for publication. 

ULF Antennas: We have extended our consideration of the excitation 

of Alfv4n waves at ULF by means of ground-based electric and magnetic 

24 

dipoles. Our earlier work showed that a field of 1 mV was produced 

at a horizontal distance of 1 km from an electric dipole of strength 

4 

X 10 coul-m, whereas with a comparable magnetic dipole, of strength 


> 


V were produced. These results sere 


2 X 10 1 ' amp-m , only 3 X 10 
based on the fields generated by waves In the earth-ionosphere wuvegulde. 
Our latest work shows that one must also consider the lateral (l.e. 
continuous eigenmode) wave for the mugnctlc dipole case, and thut this 
wave gives rise to a field strength of about 1 my . Thus, the magnetic 
dipole excites ULF waves at the earth's surface which are fully comparable 
to those generated by un electric dipole. These revised results have been 
Incorporated Into the manuscript based on Kef. J?4, which hus been submitted 
for publication. No further studies of this method of ULF excitation 
are envisaged. 


B. Nonlinear Wave Propagation 

The succors of perturbation theory in predicting small-signal 
plasma wave dispersion characteristics, has opened the way to more 
sophisticated analyses of nonlinear wave-wave and wave-particle Inter- 
actions founded on the same basic equations, l.e. Maxwell's equations 
and either a microscopic, a macroscopic, or a cold plasma description of 
the charged particle dynamics. Such analyses are frequently at the 
limits of tractablllty, and either computer simulations or improved 
analytic formalisms are required to relieve the difficulties. We have 
studied both approaches extensively under this grant, and have brought 
several significant projects to a conclusion during the reporting period: 
Compu ter Simulation : Plasma simulation on the computer entrains 

ma. , problems of Its own, in particular those of expense If a large 
number cf charged particles are to be followed, and unwanted fluctuations 

if the number of particles is small. To avoid them, we have applied a 

28 

hybrid method originally proposed by Denavlt, for effectively collislon- 
less plasmas, which provides a virtually noise-free plasma. This method 
enables us to perform simulations at signal amplitudes small enough to 
satisfy the assumptions made in the theory: noise can be reduced by six 

orders o 1 magnitude at computational costs within about a factor of two 
of previous methods. As an example, we have studied linear Landau 
damping numerically at a ratio of wave energy/plasma thermal energy four 
orders lower than in any previous work known to us. 




This work *a» followed by simulations of lurge-umpli tude longitudinal 
wove propagation. The wave distorts the time-averaged electron velocity 
distribution Into double-humped { m, which Is unstable to sideband growth 
near the phase velocity of the large amplitude wave. In practice, side- 
bands would grow from noise. In our simulations, a smull-umplltude test 
wave was Injected, and Its growth rate was measured. We lound that the 
time evolution of the distortion of the spatially averaged distribution 

function agrees with analytical results developed previously under this 

30 

grunt, and that the distortion of the averuged distribution function 
is responsible for the behavior of the sideband wave. 

In addition to the sideband growth, additional satellite frequencies 
are generated involving other sidebands and the second harmonic of the 
large-amplitude wave. A plausible four-wave interaction mechanism 
Intended to explain this was examined as part of our work under this 
grant, but j oved inadequate to account for the amplitude of the 

satellite. ^ A new mechanism was then considered, involving modulation 

29 

of the large amplitude wave. This proved highly successful in predicting 
the satellite amplitude. 

During the reporting period, a conference presentation has been 

given, ' and two manuscripts on the simulation work have been submitted 

11 12 

for publication. ’ The first of these deals with the hybrid simulation 
method, treating both linear and nonlinear Landau damping as examples. 

The second deals with sideband and satellite growth. Our development of 
simulation techniques Is now considered to be completed. We expect to 
use them on such problems as the velocity spreading of beams due to 
beam-plasma interaction. 

Lagrangian Formalism : Work has been carried out at Stanford under 

this grant, over the last four or five years, to develop and use 
appropriate Lagrangian densities for the microscopic, macroscopic, and 

cold plasma formulations. The work is described in detail in three Ph.D. 

32 -34 , 

theses, the most recent of which was completed in June 1974 and dealt 

34 

mainly with the macroscopic formulation. Although no new work on 
Lagrangians has been carried out during the reporting period, preparation 
of some of the thesis material for publication has been undertaken. In 
particular, a manuscript on the inverse Droblem of the calculus of 
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variations, i.e. derivation of an appropriate Lugrunglan from a given 
set of equations, has been written. One paper, comparing experimental 
results on nonlinear three-wave Interactions with theoretical predictions, 

derived using Lugranglan methods, has appeared du.lng the reporting 

7 

period, and a conference presentation has been made on the linear 

4 

problem of resonances In an Inhomogeneous plasma, l.e. Tonks-Dattner 
resonances . 

C. Long belayed Echoes 

From 1967 to 1973 , nn Ionospheric sounding program was carried out at 
Stanford to determine whether u curious phenomenon reported by van der Pol 

ts 

and Stormer in 1928 really occurs or not, and if so to determine Its 
origin. The effect is slip y that radio signals, e.g. Morse dashes can 
occasionally return from the Ionosphere with tens of seconds delay, 
rather than the few ms delay expected. Sporadic observations by a 
large number of amateurs lend credence to the Idea that the effect is 
genuine, but no systematic studies have been unequivocally successful 
In demonstrating its existence. Some suggestive echo data from the 
Stanford experiments have U'd to a theory of the phenomenon based on 
beam-plasma Interaction due to precipitating fast particles Interacting 
with the Ionosphere so as to allow low group velocity propagation (hence 

35 

long delay) without undue colllslonal attenuation. 

Last year, NSF support for the LDE program came 1'* an end, and was 
replaced by support from this NASA grant. The observational program 
has not been continued, but we have strong reasons for pursuing the 

theoretical work: first, our most recent numerical results show very 

2,8 

encouraging agreement with the observed characteristics of LDE, and 
second, the Space Shuttle offers the possibility of studying the 
phenomenon from the topside of the ionosphere. On the topside, there is 
less radio noise to contend with, and frequency usage is not constrained 
by FCC allocation. Consequently, much more comprehensive measurements 
should be obtainable rupldly. 

During the reporting period, a Ph.D. thesis has been completed by 

8 

Sears reviewing the results obtained at Stanford and elsewhere. The 
thesis also contains some preliminary steps towards an analysis of the 


phenomenon. It 1 b postulated that ordinary mode energy from the 
transmitter couples linearly Into a longitudinal plasma wave, propagating 
with low group velocity («*» 1 km/a) along the earth's magnetic field lines, 
close to the ordinary mode reflection height. If there are some hlgh- 
er.ergy («*- keV) electrons precipitating along the field lines, they may, 
through the mecnunlsm of beam-plasma Interaction, offset the considerable 
colllslonal damping that a low group velocity wave would otherwise suffer 
over times of the order of 10 s. The plasma wave may then couple linearly, 
for example, at a local plasma lnhomogenelty, to the ordinary mode, which 
can propngate downwards to the receiver. 

The Initial calculations by Sears for plausible assumed high energy 
electron fluxes have two defects. First, the maximum Integrated group 
delay Is about 1 s, l.e. i.n order of magnitude too low, and second, the 
growth that the signal wave-packet must experience over some parts of 
Its path, to compensate for damping over others, would drive the Inter- 
action into nonlincur saturation. Thus, a modified mechanism capable of 
predicting weaker growth and longer delay must be sought. During the 
reporting period, we have considered two additional factors: Ionospheric 

drift, and plasma wave propagation oblique to the earth's magnetic field. 

It. was hoped that the effect of drift might be to Increase the group 
delay appreciably, and possible to allow coupling out of the Ionosphere 
to occur at the reflection height of the extraordinary mode. It was 

found, however, that plasma drifts 1 km/s would be required, l.e. 

36-jb 

appreciably larger than those reported in the literature. Our 

studies of growth rates and group delay for oblique propagation have 
only begun recently. It is anticipated that they will be concluded 
during the coming reporting period. 

I). Pulse Propagation 

Although the propagation of continuous waves through essentially 
homogeneous plasmas has been studied for over forty years, it Is only 
in the last decade or so that transient propagation has received much 
attention. Under this grant, we have already conducted some pulsed 

39 

propagation studies of cyclotron harmonic waves, and obtained suffi- 
ciently good agreement with the theoretical dispersion characteristics 
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for the group delay to be used as a non -perturb Inn plasma diagnostic 
technique. With the Space Shuttle, the opportunity Is offered of cutrv.ng 
out such experiments under conditions whjre Instrumentation problems are 
very much simpler to solve than In the laboratory. It Is our Intention 
to examine the range of possible experiments. To begin with cold plasma 
propagation, packets In the ordinary and extraordinary modes could be 

excited with effectively delta-function pulses (*• 10 ns), and exhibit 

4 0 

the features predicted by lirillouln and Sommerfold around WW I. These 
still await detailed verification with respect to such phenomena » 
"forerunners". Warm plasma wave propagation in the Landau and cyclotron 
harmonic modes could also be studied. 

So far, our effort has been devoted to analyzing the wave-packets 
produced by del tu-f unction excitation of a cold magnetoplasma. This 
shou’ v- ' vo two branches of right-hand polarized propagation, and one 
lrf nand polarized branch, for wavevectors parallel to the earth's 
magnetic field. The analytical results must be studied numerically to 
make their significance apparent. This has been done duiing the reporting 
period. The numerical results allow us to identify the right- and left- 
hand polarized branches of the dispersion relation. For these waves In 
a weak magnetic field ( < 0.1), an approximate analytic solution 

has been found. 

The whistler branch of the right-hand polarized mode cun be 
distinguished In the numerical solutions. To first order, the velocity 

of propagation for this branch corresponds to the classical group 

1/2 

velocity, v^/c ~ ( 27 /' *0 (‘ c / u ) • The detailed structure of the 

forerunners to this wave Is effectively a rapidly Increasing sinusoidal 
wave that merges into the main wave-packet. An analytic expression for 

the arrival of the whistler wave-packet Is currently being developed, 

4l_4 4 

using an asymptotic method, ~ J and should provide more detailed 
Information on the structure of the forerunners. Preliminary results 
indicate that, for typical ionospheric parameters at heights of about 
k^O km, this expression for the whistler branch will be accurate to 
within 1^ when the transmitter and receiver are separated by a distance 
of a few hundred meters or more. 


9 
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E. Ionospheric Heating and Backscatter 

Over the last few years Ionospheric heating experiments have been 
carried out with the ESSA high-power (<*» 1 MW) sounder, near Boulder, and 
at Areclbo (Ref, I gives a brief review of the literature). Under such 
conditions, l.e. when a hlgh-lntenslty signal propagates through the 
Ionosphere, three separate nonlinear mechanisms can operate to Increase 
the "resistive" losses above those due simply to the colllslonal effects 
included in linear theory. The mechanisms are: (a) thermal 

Instability, (b) oscillating two-stream Instability, and (c) parametric 
Instability, The coupled-mode formalism that we have used In our wave- 
wave analyses ’.'nder this grant provides a suitable approach to the last 
two mechanisms. With its aid, both of them can be shown to depend on 
the nonlinear Interaction between ordinary, Langmuir, and lon-acoustlc 

1*9,50 

waves. 

Such instabilities can be, and have been, studied by the radar 
bacXscatter technique. This .aises an important question: though the 

heating signal may be effectively monochromatic, Langmuir and Ion acoustic 
waves will be excited over a band of frequencies; what will be the shape 
of the backscattered spectrum? Solution of this problem involves 
taking ac air.t of the nonlinear saturation mechanisms which limit the 
parametric growth. During the reporting period, we have removed some 

of the shortcomings of a highly simplified approach developed under 

51 

this grant about two years ago, and have obtained satisfactory predictions 
of the locations and separations of the main lines in the backs atter 
spectrum under conditions appropriate to Areclbo. This analyses has 
been written up recently for publication,^ and is regarded as completed. 
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A. Whistlers 

We wish to pursue the project on the excitation of VLK waves by 
charged particle beams described in Section II A. In particular, it 
would be desirable to calculate the total power radiated by the beam, 
rather than Just that part which is ducted. This would provide Important 
information for experiments where receivers are to be located in tne 
vicinity of the Shuttle. This calculation is not possible within the 
cold plasma analysis that we have used up to the present, since the 
fields become infinite at a certain angle, giving rise to an infinite 
radiation resistance. This is the familiar problem in antenna theory of 
the "resonance cone", and always occurs for waves which have open wave- 
normal surfaces. The difficulty can be removed by introducing electron 
temperature into the model, thereby giving rise to convergent integrals. 
We shall extend out analysis of the problem, by using warm plasma theory, 
in order to solve this important problem. 

So far, we have considered a circular array in which all of the gunr 
are mount d on the circumference and fired sequentially, bowden's model 

distributes the guns over the circular area, and all are pulsed 

25 

simultaneously. We shall extend our theory to cover this case. It 
is expected that the power radiated will then be increased by about 
one order of magnitude. 

U. Nonlinear Wave Propagation 

It is anticipated that two papers will be written during the next 

six months on the Lagrangian formalism for the macroscopic plasma 

34 

description developed in the thesis by Peng. The first will describe 
the derivation of the appropriate Lagrangian and Hamiltonian densities; 
the second will consider perturbation expansions of the Lagrangian to 
obtain linear wave propagation and nonlinear wave-wave- coupling coeffi- 
cients. No further work on the formalism itself is envisaged: it will 

simply be used as occasion demands in our analyses ol plasma wave and 
resonance phenomena. 
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C. Long Delayed Echoes 

Farther work is required to test the beam-plasma mechanism outlined 
in Section II C, and described In detail in i<ef. 8 . In particular, 
computations of group delays and growth rates for propagation oblique 
to the earth's magnetic field lines ere required. A major publication 
will then be prepared embodying these results and those contained in 
kef. 8. 

D. Pulse Propagation 

When the pulse propagation phenomena have been thoroughly studied 
for propagation parallel to the earth's magnetic field, including effects 
of positive ions and collisions, the analysis will be extended successively 
to the limiting case of perpendicular propagation, and the general case 
of oblique propagation. It will be of particular interest to determine 
which features of the wave-packet can be observed as measures of the 
ionospheric parameters (density and magnetic field). 

E. Ionospheric Heating and Backscatter 

In Section II E, we noted that ionospheric heating experiments had 
been carried out using ground-based transmitters located near Boulder 
and at Arecibo, and had yielded evidence of the excitation of nonlinear 
wave interactions. The Space Shuttle would alljw a more modest trans- 
mitter antenna system to be carried, but still capable of providing 
similai power fluxes in the ionosphere. Since radar backscatter 
represents a powerful technique for studying the nonlinear interactions 
likely to result, involving transverse ordinary and extraordinary waves 
with longitudinal electron plasma and ion acoustic waves, we shall spend 
some time during the next reporting period investigating the feasibility 
of carrying an adequate backscatter system on the Shuttle to make both 
coherent and incoherent backscatter measurements. 
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